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Abstract 

Solutions of the scalar Helmholtz wave equation are derived for the analysis 
of the transport and thermodynamic properties of the two-dimensional disk 
and three-dimensional infinitely long straight wire in the external uniform 
longitudinal magnetic field B under the assumption that the Robin bound- 
ary condition contains extrapolation length A with nonzero imaginary part 
Aj. As a result of this complexity, the self-adjointness of the Hamiltonian is 
lost, its eigenvalues E become complex too and the discrete bound states of 
the disk characteristic for the real A turn into the corresponding quasibound 
states with their lifetime defined by the eigenenergies imaginary parts Ei. 
Accordingly, the longitudinal flux undergoes an alteration as it flows along 
the wire with its attenuation/amplification being independent too. It is 
shown that, for zero magnetic field, the component Ei as a function of the 
Robin imaginary part exhibits a pronounced sharp extremum with its mag- 
nitude being the largest for the zero real part A r of the extrapolation length. 
Increasing magnitude of A r quenches the Ei — Aj resonance and at very large 
A r the eigenenergies E approach the asymptotic real values independent of 
Aj. The extremum is also wiped out by the magnetic field when, for the large 
B, the energies tend to the Landau levels. Mathematical and physical inter- 
pretations of the obtained results are provided; in particular, it is shown that 
the finite lifetime of the disk quasibound states stems from the Aj-induced 
currents flowing through the sample boundary. Possible experimental tests 
of the calculated effect are discussed; namely, it is argued that it can be 
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observed in superconductors by applying to them the external electric field 
£ normal to the surface. 

Keywords: Robin boundary condition, quantum disk, Ginzburg-Landau 
theory of superconductivity, wave scattering, quasi bound states 



1. Introduction 

A solution \I/(r) of the scalar Helmholtz wave equation 

V 2 ^(r) + fc 2 tf(r) = (1) 

is strongly influenced, for the spatially confined domain i?, by the require- 
ments imposed at the sample interface dQ. A demand zeroing at the edge a 
linear superposition of both \l/(r) and its spatial derivative is known as Robin 
boundary condition |l( with its most general form written as 



(2) 

dQ 



n being an inward unit vector normal to the confining interface dQ and 
Robin parameter 1/A being, in general, a function of the one- (ID), two- 
(2D) or three-dimensional (3D) radius- vector, A = A(r). Such boundary 
conditions with real A emerge naturally in different physical and chemical 
systems wha t sp urs their theoretical investigation by mathematicians and 
physicists [2 46 . Examples of the vibrating strings and membranes are so 



well known that they are discussed in the undergraduate textbooks [3]. A 
time-independent Green function method involving a multiple reflection ex- 
pansion was used for calculating the distribution of eigenvalues of ([I]) for the 
arbitrary domain with the Robin condition on its surface |4|. Extension to 
the vector fields and to an arbitrary number of dimensions followed soon j^] 
(errata of these two papers are fixed in The requirements of Eq. (J2J) 

are used in thermodynamics and statistical physics for the description of the 
properties of the ideal gases enclosed in the multiply connected bounded con- 
tainers |l9| . Static and dynamical Casimir effects [47j were calculated on the 



straight plates and miscellaneously curved surfaces |1(J, ll5|, l23|, |3JJ, |42|, |43 
supporting the Robin condition, Eq. fl2]). It was also used in the study of 
the correspondence between a field theory on anti-de Sitter space and a con- 



formal field theory on its boundary [13] . In solid state physics, for the de 



2 



Broglie particle with mass m p the wave vector k in (pQ) is expressed as 

k = ^/2m p E/h, (3) 

where E is a total energy of the particle and H is the reduced Planck constant. 
The electron in semiconductor nanostructures is satisfactorily well described 
by the zero Robin distance, A = 0, what is a limiting Dirichlet case of the 
general demand of Eq. (T5]) , 

*I M = 0- (4) 
On the other hand, in the phenomenological Ginzburg-Landau (GL) the- 



ory of superconductivity [48|, |49[ the Cooper pair behavior is described by 
the order parameter \l/(r) with its square being equal to the density of the 
superconducting particles n s , 

^ = |^(r)| 2 . (5) 
In general, the order parameter satisfies nonlinear GL equation 

V 2 ^(r) + A; 2 ^(r) + |^(r)| 2 *(r) = 0. (6) 
Here, the role of the energy E in ([3]) is played by the parameter —a: 

with T being the actual temperature of the superconducting material, T c 
being the bulk critical temperature at zero magnetic field, £ (0) representing 
zero-temperature coherence length, and Cooper pair mass being equal to the 
double bare electron mass m e . An expression for the second GL parameter 
(3 in Eq. (E]) is derived from the microscopic equations of superconductivity, 
and for the "pure" superconductors it reads 

l h 2 l 

^iV(O) 2m p £ 4 (0) {kBTc f {} 

with (0) being the density of states at the Fermi energy and ks the Planck 
constant. For the 'dirty' superconductors (alloys) £ 2 (0) in GL parameters a 
and (3 should be substituted by £ (0) /, with / being the mean free path of 



the material [49|. In the temperature range close to the transition to the 



normal state, the order parameter is small, and the cubic term can be safely 
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neglected 49[ thus simplifying the nonlinear GL equation into its linear 
Helmholtz counterpart, Eq. (pQ), which is much easier to handle. As it follows 
from OH]), the influence of the cubic term can be also strongly suppressed 
by the appropriate change of the materials 45|, |50| which leads again to the 
analysis of (JT]) whose smallest eigenvalue E min , according to ([7]), defines the 



the sample temperature [49 



T = T, 



1 - 



2m p e (0) 

h 2 



E„ 



(9) 



This expression shows an enormous importance of minimizing the lowest 
eigenenergy E min since this leads to the higher superconductor tempera- 
tures. One of the possible methods of this minimization is the curving of the 



previously straight sample [50H53 



In the GL theory the order parameter behavior near the surface is gov- 
erned by the same Eq. (j2J). The de Gennes distance, as is the length A called 
in this field of physical research, is a quantitative measure of the interaction 
between the superconductor and the ambient environment. For example, the 
border with the vacuum or the insulator is described by the infinite extrap- 
olation length, 1/A = 0, which leads to the Neumann limit of Eq. (j2j), 

It = °- (10) 

For the superconductor border with the normal metal the de Gennes dis- 
tance takes the values in the interval < A < oo where the limit of Eq. (j3J) 
is approached for the bordering ferromagnets. Solutions of ([6]) were found 
and analysed for the different domains Q and Robin factors A 54M571]; in 
particular, a crucial role of the extrapolation length in the nucleation of the 
superconductivity in thin films was underlined. Interestingly, the de Gennes 
distance can take negative values, A < 0, for the boundary with the super- 
conductor with higher critical temperature T c . This theoretically predicted 
enhancement |5ll.l52l.l58| was indeed observed in cold worked Ino.993Bio.007 foils 
581 ] and tin samples 59(. It was estimated that in the former case the ex- 
trapolation length is A ~ — 1 /im [5l| . Negative extrapolation length appears 
also in the theory of twinning plane superconductivity 60] . Robin boundary 
conditions, Eq. (j2J), are also utilized in the study of magnetic multilayers 45 



and ferrite- filled waveguides [6ll465l|. Recently, they were also used to model 
the tunneling current and the optical gain of GaAs/ALGai^As quantum 
cascade lasers 44 . 
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To take into account the influence on the charged carriers of the external 
magnetic field B, one needs in a standard way to incorporate it into the 
Helmholtz equation through the introduction of the vector potential A, 



B = V x A. 



Then, for the de Broglie particle, Eq. ([I]) turns to 



V -4A)~^(r) + A; 2 #(r) 
n 



0. 



(12) 



where q is the particle charge; in particular, for the Cooper pair, q = — 2e, 
with e being the absolute value of the electronic charge. Eq. (jH]) is modified 
in the same way: 



(V - ^A) 2 *(r) + fc 2 *(r) + |^(r)| 2 (r) = 0. 
In the GL theory, the boundary condition changes to 



n 













dQ A 



(13) 



(14) 



Theoretical analysis of the combined influence of the magnetic field and real 
Robin length on the superconductors [6614891] revealed, among other features, 
a significant influence of the parameter A on the nucleation of the supercon- 
ductivity, critical magnetic fields and localization properties of the supercon- 
ducting state. 

A natural extension of the real Robin distance is its continuation into the 
whole complex plane A = (A r ,Aj). For example, the boundary conditions 
with purely imaginary A = zAj for the ID and 2D geometries are extensively 
used 9CH94 | in the actively developing field of PT-symmetric quantum me- 



chanics 95l-l97j. Apart from the purely theoretical interest [98|, the analysis 
of the systems with the complex Robin lengths has large practical applica- 
tions as they model scattering phenomena in different media j99j. Consider, 
for example, electromagnetic or acoustical processes when the wave vector 
k from ([!]) is the ratio of the actual frequency u and the speed c of the 
propagation of the corresponding oscillations: 



k 



UJ 

c 



(15) 
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Figure 1: Schematic of the infinitely long straight superconductor wire of the radius po 
subjected to the uniform magnetic field B parallel to its axis. Channel walls support 
Robin boundary condition, Eq. ([2J, with complex extrapolation length A uniform along 
the wire surface. Cartesian (x, y, z) and cylindrical (p, (p, z) systems of coordinates are also 
shown with their origins lying on the waveguide axis. Curved arrow shows the azimuthal 
direction in which the polar angle ip grows. 



It is well known that a porous lining of the walls of the soun d duc t results 



in nonzero and, in general, complex acoustical admittance Y 100] . Corre- 
sponding boundary condition for the air pressure p that satisfies Eq. <^ is 
then describ ed by (El) with the coefficient 1/A being prop ortio nal to the ad- 



mittance Y [10lHl09l | . The same boundary requirements [110] appear in the 
study of the propagation and scattering processes in the impedance electro- 
magnetic waveguides with the fields independent of the tange ntial direction 



1111. In this 



when the division onto the TM- and TE'-modes is preserved 
case, the coefficient 1/A in (j2J) is either proportional or inverse proportional 
to the complex surface impedance w an d its imaginary part is determined 
by the real component of the impedance 111] . Apparently, the same model 
with complex extrapolation length A describes the processes in the ferrite- 



filled guiding structures with complex permeability [6l H63l Ill2] . Dissipative 
Schrodinger operators with complex ID Robin boundary conditions were 
analyzed for the study of the drift-diffusion phenomena in low-dimensional 



semiconductors [113H116] . Superconducting de Gennes distance can take 
complex values too. Namely, it was argued recently 1171 . Ill8] that the elec- 
tric field £ applied perpendicularly to the surface should be accounted for 
in the total de Gennes distance A to t by the addition to the inverse zero-field 
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extrapolation length 1/A of the extra term S/U, 

1 1 



A to t A 



+ £ (16) 



with the potential U s being expressed through the parameters of the GL 
theory: 

1 _ w 2 d ^T c \q\e s 
U s d In n s rripC 2 

Here, the reduction factor r\ is the ratio of the band gap extrapolated to the 
surface to the gap value at the surface; dimensionless GL parameter k is the 
ratio of the zero-temperature London penetration length A(0) to the coher- 
ence length, k = A(0)/£(0); e s is superconductor ionic background permit- 
tivity, and c is speed of light. According to the estimates [49j], the coefficient 
1] is of the order of unity. E lectr ic field influence on superconductors is still 



an open question, see Ref. |119| and references therein. If the permittivity 
e s has a noticeable imaginary part, so does, according to f|T6|) and ffTT|) . the 
total extrapolation length too. It is well known that the real value of the 
de Gennes distance is due to the fact that no superconducting current flows 
through its surface |49j : indeed, as we shall show below, a complex part of 
the extrapolation length forces the current to acquire nonzero normal (with 
respect to the interface) component. 

Thus, solutions of the Helmholtz equation with complex Robin parameter 
are applicable in many areas of physics. Even though some investigations 
have tackled different aspects of this problem, as the above review demon- 
strates, so far no comprehensive detailed dependence of the bounded domain 
properties on the extrapolation length A was given and no magnetic field 
influence was addressed. Present research closes this gap for the simplest 
case of the 3D straight channel and its cross-sectional 2D disk; namely, the 
problem of the infinitely long straight waveguide of the circular shape with 
the boundary conditions, Eq. fl2]), supporting complex extrapolation length 
A is formulated and solved. Uniform magnetic field B is pointed in the wire 
axis direction. As stated above, without magnetic field such a model and its 
analysis based on the solutions of Eq. ([[]) is relevant for the description of 
the wave dynamics in lined acoustical and impedance electromagnetic ducts 
as well as ferrite or superconductor wires. Turning on the field confines the 
applicability of the model to the latter situation when Eq. ffl2|) comes under 
scrutiny. In either case, analytical solutions of the corresponding equation 
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are derived and analyzed. Based on these solutions, currents and magneti- 
zations are calculated in the whole complex de Gennes plane A r — Aj. It is 
shown that the magnetization for the purely imaginary extrapolation lengths 
exhibits basically the same dependence as for the real Robin distances while 
the longitudinal current along the wire shows resonance features in its de- 
pendence on the imaginary A. The physical and mathematical reasons for 
such highly nonmonotonic behavior are described in terms of the transverse 
currents induced by the imaginary component of the extrapolation length. 
These currents flow through the superconductor boundary and in this way 
they change the number of the Cooper pairs participating in the longitudinal 
transport. Possible experimental confirmations of the developed theory are 
discussed. 

The paper is organized as follows. In Section II our model is presented 
and a necessary formulation of our method is given. Section III is devoted to 
the presentation and detailed mathematical and physical interpretation of the 
calculated results. Summary of the research with possible future extensions 
is provided in Section IV. 



2. Model and formulation 

Infinitely long straight wire of the circular cross section with the radius 
Po is placed into the uniform magnetic field B with its direction coinciding 
with the channel axis (Fig. [T]). Cylindrical walls of the waveguide support 
boundary condition, as described by Eq. ([2]) with the uniform along the 
length and circumference extrapolation length A. We do not confine the 
value of A to be real concentrating on the properties of the structure at the 
complex de Gennes distances. Geometry of the system dictates a natural 
choice of the cylindrical system of coordinates (p, <p, z) with its origin lying 
at the circle center and the z axis being parallel to the waveguide. We will 
seek the solutions of Eq. f[T2"l) with the vector potential from Eq. (fTTj) written 
in the symmetrical gauge, 

A = (0, Bp/2,0). (18) 

The fact that we treat Eq. (fT2|) with the uniform field means that, for the case 
of superconductors, we restrict our consideration to the range of the magnetic 
intensities and temperatures close to the transition between superconducting 



and normal states [49j even though the results obtained are covered by the full 
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Figure 3: (a) Real Enm and (b) negative double imaginary T nm parts of the transverse 
energy E^ m as functions of for A r = and zero magnetic field for several radial n and 
azimuthal m quantum numbers. Solid lines depict the state with n = m — 0, dotted lines 
are for n = 1 and m = 0, dashed lines - for n — 2 and m = 0, dash-dotted lines - for 
n = and |m| = 1, and dash-dot-dotted lines - for n — and |m| = 2. The inset in panel 
(a) shows real energy Ego as a function of the real extrapolation length A r (dotted curve) 
superimposed on the dependence of the real part Eq of the complex energy Eq versus 
purely imaginary extrapolation length iA,; (solid line). 
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GL theory (see below) and have much wider validity range. We will operate 
with the energy E through which the wave vector k is expressed, according to 
(j3J). Such a treatment describes a superconductor wire. A transition to the 
frequency u of the acoustical or electromagnetic oscillations can be readily 
done with the help of (1151) . It is important to remark here that, for the case 
of ferrites, the static magnetic field influence on the propagation properties 
is through the dependence of the magnetic permeability tensor /ion5 and, 
thus, one needs to consider Eq. (pQ) where the speed of propagation depends 



on the magnetic field through /i [64], |65j. In order to make the obtained 



results as generic as possible, it is convenient to express all quantities in 
dimensionless units. In the GL theory, one usually chooses as a unit of 
distance the coherence length £(0); however, to make our results directly 
applicable in the acoustics and electrodynamics, we choose the radius po as 
a natural unit of length; accordingly, if not stated otherwise, all energies will 
be measured in units of the ground-state energy ir 2 h 2 j{2m v pfy of the infinite 
Dirichlet ID quantum well of width p ; all momenta, in units 1/poi magnetic 
fields, in units of h/(\q\p%); magnetization, in units of h\q\/(2m p ); 2D current 
density, in units of gft/(m p pg); current, in units of g^/(m p p ); time, in units 
of 2m v p\j{j: 2 fi). In these units Eq. ([3]) transforms to k = iryE. Then, 
Eq. (fl2~|) for the order parameter \I/(p, <p, z) becomes 

1 d f d^\ , I d 2 ^ : 1 2 2t , , d 2 ^ 



pdp \^ dp ) p 2 dip 2 * dip 4 ^ dz 2 0(9) 
with E being a total energy of the particle. Factoring out the z-dependence 

*(p,<p,z) = e ik "1>(p,<p) (20) 
leads to the equation for the transverse function ip{p, <p): 

1 d ( _d$\ , 1 d 2 ifj , 1 2 2 2 ± 



- P T P [»t ) + lB t - 4 B " v + * E i ' = °- (21) 

Longitudinal wave vector fc z and the transverse energy i? -1 are related as 

k z = nVE - E x . (22) 

Rotational symmetry of the system allows one to separate out the angular 
and radial dependencies in the transverse function ip(p, <p): 

?/w(p, <p) = ^ e im ^„ m (p). (23) 

V 27T 
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Here, m = 0, ±1, ±2, . . . and n = 0, 1, . . . are the azimuthal and the principal 
quantum numbers, respectively. In this way one arrives at the equation for 
the radial function R nm (p): 

d 2 R nm 1 (IRnrr, ( TTl 1„ V , Jpl p _fl (OA\ 

tinm + 7T £j nm K nm — U. (24 ) 



dp 2 



I 1 dR nm 
p dp 



m 1 

7 + 2 B " 



This equation is supplemented by the boundary condition for the function 
Rnm(p)- Our choice of the vector potential in the form of ffl"8]) drops it out 
from (fl"4"l) which becomes 



dp 



0. 



(25) 



Eqs. f[2"4"|) and f[2"5l) constitute the problem of finding the eigenfunctions 
Rnm(p) and eigenenergies E^ m of the 2D circular disk with its circumference 
supporting the boundary condition with, ingeneral, complex A. We mention 
that the field- free case of the real constant [i, 37] or varying 37, 39, 4^] along 
the perimeter extrapolation length was considered before, as was the case of 



the quantum dot with real constant A in the uniform magnetic field [69h71 
with the Dirichlet 1 1 2 0M 1 2 5f l or Neumann [71 

(E 



126-1301 limits. 



Analytical solution to Eq. 
hypergeometric function M(a, b; x 



is ex pressed via the Kummer confluent 



131 



Rr, 



xM 



(p) = Inm exp 

m + \m\ + 1 



V 

4 H 



V 

2 H 



\m\/2 



7[- 



E 



\m\ + l--Bp< 



(26) 



2 IB 

with the real coefficient 7„ m determined from the normalization condition 
which is either of the form from (JHJ) for superconductors or 

-i 



\Rnm{p)\ pdp 



(27) 



for the electromagnetic and acoustical waveguides. Applying boundary re- 
quirement, Eq. (I25p . one arrives at the following transcendental equation for 
the determination of the energies E^ m : 



-BM' 



B 1\ 
- + - M 
2 k) 

m + \m\ + 1 



m + I ml + 1 



7T 



Ei 



71 



2 

2 TP- 



2 B 



\m\ 



B 
1; 2 



2 B 







\m\ - 





(28) 
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with prime denoting a derivative of the function with respect to the last 
argument. Instead of the derivative, one ca n use the Kummer function itself 



with the different parameters according to [131 



M'(a, b, z) = -M(a + 1,6 + 1,*). 
b 



Utilizing properties of the Kummer function [1311 ]. it is elementary to show 
that in the limit of strong magnetic field, B — > oo, solutions of (|28|) tend 
from below to 

2 ( m + |m|+l\ „ 



B 



which are familiar Landau levels jl32 , 133 ]. In the opposite limit of the 
vanishing magnetic field, B — > 0, Eq. f)28p transforms to 



^\^Li J \m\ (^^EL) + J J \m\ (^V^m) = ( 30 ) 

This equation was 
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with J m (x) being a Bessel function of the order m 

derived before for the case of the real @, E 0, H HEl and complex [l0^ 
extrapolation lengths. Zero-field wave function is expressed then as 



Rnm(p) = lnmJ\m\ ^7Ta/ 



(31) 



Of course, this last expression could be obtained directly from ( 12^1) if one 
puts B = in it and, subsequently, Eq. ( 130]) is derived from it by means of 
( 125]) . From (128]) and ( 130]) it straightforwardly follows that for the complex 
extrapolation lengths, A = A r + iAj with real A r = Re(A) and Aj = Im(A), 
the energies are complex too, 



E~ 



with real E n r m and T nm , and the following property holds: 



E± m (A) = E± m (A) 



(32) 



(33) 



with the overline denoting a complex conjugate value. The same is true for 
the corresponding radial functions too. Minus sign and the fraction in the 
imaginary part of the right-hand side of (13"2"]) were taken to be consistent with 
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the scattering theory [133l |134| which teaches that an imposing of the com- 
plex boundary conditions eliminates the self-adjointness of the corresponding 
Hamiltonian what results in the complexity of its eigenenergies. Depending 
on the type of the requirement at the border (in our particular case, the sign 
of Aj), the imaginary part can be either positive, T < 0, or negative, T > 0. 
In the latter case, the true bound level with the discrete at Aj = energy 
transforms into the quasibound state broadened over the width Y around its 
central value E^'\ The particle does not stay forever in such quasibound 
state but after the time ^ 

J- nm 

leaves the system. As ( )34l) shows, the lifetime r nm is inversely proportional 
to the half width. The negative value of V means, in turn, a flux of the 
particles into the system and their accommodation inside it with the temporal 
accommodation rate being proportional to l/|r|. We will discuss these two 
situations in more detail later in the text. 

It is important to state that even though f[2"8"j) and f[3"Uj) were derived for the 
linearized GL equation, it can be shown that for the thin superconducting 
disks they follow also from the complete nonlinear GL theory. For doing 
this, one needs to employ the method of minimizing of the free energy F of 
the superconducting state with respect to the coefficient linking the order 



parameter of the nonlinear equation with its linear counterpart, Eq. (|23J) . As 
this procedure is described in detail for the superconducting rings 



81 



83 



we do not repeat it here referring the reader to the original research, Refs. 



81, 83 



Knowledge of the eigenfunctions and eigenenergies allows one to find other 
physical characteristics of the structure. For example, in the case of super- 
conductors a magnetization operator 

M = i|_I V (35) 

is used to calculate the magnetic moment M z of the 2D disk: 

M z = (if) nm (p,(p)\M\if) nm (p,ip)). (36) 

One immediately gets: 

\R nm (p)\ 2 P 3 dp . (37) 





b r 1 








~2Jo 
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It is worthwhile to note that the zero -field case allows analytical calcula- 



tion of the integrals in (I27p and (1371) [135l . Il36| for the real extrapolation 
lengths when the radial f unctions are real too; however, since there are no 
in known to us literature 13 lL ll35Ml38| analytical expressions for these inte- 
grals with the complex functions R n m{p), their direct numerical quadrature 
was performed. The same procedure was used for the integrals involving the 
Kummer functions for the nonzero magnetic fields with either real or complex 
extrapolation length. 
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An expression for the superconductor current density j is written as [48 . 



j = Im ( r ) (r)] + (r) * (r) 



(38) 



Similar formula (wit h, of course, A = 0) can be used for the Poynting vector 



in el ectrodynamics [1391 ] or for the sound energy density flux in acoustics 



140| . In our configuration, for the 2D quantum dot without the z dependence 



it transforms to 



1 

2^ 



Im 



-dR\ 



m 1 
- + -Bp 



I D|2 

\R\ e u 



(39) 



e p and e v along the radial and azimuthal 



with the unit orthogonal vectors 
directions, respectively. The real value of the extrapolation length A guar- 
antees that no supercurrent flows through the border [49|. This general 
statement can be checked directly with the help of (|38|) and (JI4"|) which in 
our dimensionless units reads 



n(V + iA) * 



an 



1 

A 



(40) 



an 



For our system, Eq. (|39|) vividly proves it since the normal component of the 
current described by the term at the radial unit vector, for the real radial 
functions is an identical zero at the cylindrical surface, what is the case for 
the real Robin coefficient, as Eq. ( 1251) demonstrates. 

For the wire the total longitudinal supercurrent J z is given as 



J z 



2tt 



pdp j dipj z = Re(k z ) exp [— 21m(k z )z] . 



(41) 



Eq. (jSJ) shows that the longitudinal dependence of the density n s is of the 
same form: 

n s = \R (p) | 2 exp [-2hn(k x )z] . (42) 
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Expressions for Ke(k z ) and lm(k z ) follow straightforwardly from ( 12 2 p and 

mi- 



Re(k z ) 
lm(k z ) 




E - E (r) 



+ (r nm /2) 2 + 
r 11 

x nm/ f 



E - E {r) 




(43a) 
(43b) 



l-i J^nm 



+ (r nm /2) 2 + 



rp _ rp(r) 



Note that the real and imaginary parts of the wave vector k z are not inde- 
pendent but are related to each other as 



R,e(k z )lm(k z 



7T 



r. 



(44) 



As is seen from (|4ip and (143]) . for the real extrapolation lengths when, as 
stated above, T nm = 0, the supercurrent, or the electromagnetic or acoustic 
energy flows along the channel without changes in its intensity, R,e(k z ) = 

ty\J E — En m} lm(k z ) = 0] however, as soon as the Robin distance acquires 
a complex component, the nonzero imaginary part of the wave vector forces 
the current to change as it flows down the waveguide. For the supercon- 
ducting, electromagnetic or sound transport phenomena a positive sign of T 
determines an attenuation of the corresponding energy flux and the parti- 
cle density with the fading intensity being proportional to its magnitude, as 
Eqs. (jU]) - (jl5|) demonstrate. However, in our case, as Eq. ( 133|) manifests, the 
value of T can be negative too since the sign of the imaginary part of the 
total extrapolation length A tof in our model can be changed by the simple 
switching of the gate voltage, as it follows from (Tl6|) and ( |T7I) . In this case, 
instead of dissipating the energy into the surrounding environment, as is the 
case for T > 0, the channel absorbs it from outside. In reality, the increase 
of the current approaches some saturation value since an infinite pumping 
of the energy is physically senseless. This saturation is not covered by our 
model. Moreover, in the superconductor wire the current upon reaching some 
critical value destroys the superconductivity. This, too, should be covered 
by more elaborated theories. 

Used in the present research GL theory, being phenomenological in its 
nature, does not provide an explanation of the microscopic mechanism of the 
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carrier transformation in the superconducting and nonsuperconducting me- 



dia. Invoking the Bardeen-Cooper-Schrieffer theory |49] and our earlier dis- 
cussion about the electric field influence on the Robin parameter [see ffl6|) and 
(fI7|) ]. one can say that single electrons pushed by the appropriately directed 
perpendicular electric field from, say, normal metal into the superconductor, 
are bound inside it into the Cooper pair and, thus, make an additional con- 
tribution to the total longitudinal supercurrent as the nonsuperconducting 
flow of the normal electrons is an identical zero in the thermal equilibrium. 
Switching the direction of the electric field forces the Cooper pair to move 
out of the superconductor into the medium that can not support this bound 
state with the corresponding decrease of the longitudinal supercurrent. In 
the same way, the electrodynamics operating with the complex impedance w 



llOl llll] does not provide the detailed microscopic description of the wave 
attenuation which is explained by the electromagnetic energy transforma- 
tion into the joule motion in the waveguide metallic walls. These nuances 
are simply incorporated into the value of w. 

To get more insight into the physical meaning of the complex energy, it is 
instructive to calculate the 2D divergence of the supercurrent in the quantum 
dot. Eq. ( )39|) shows that only the radial component contributes. Utilizing 
properties of the Bessel functions, one gets: 



7T , 



divj = -|i?|T. (45) 

For the waveguide this expression in each cross section z = const should 
be multiplied by the familiar from above exponential factor exp [ — 2Im(A;^) 2r] . 
Thus, the imaginary part of the energy defines the divergence of the current 
with its positive (negative) value meaning that the corresponding spatial 
point serves as a sink (source). This explains the longitudinal current change 
for the complex energies as for the positive (negative) V the Cooper pairs leave 
(enter) the superconductor through the cylindrical surface thus decreasing 
(increasing) in the wire the number of charge carriers participating in the 
transport along the z axis. 

In conclusion of this section, we note that the complex energy E means, 
according to (Q, that the temperature T becomes complex too. According 
to our earlier interpretation, real part of the temperature means an actual 
temperature of the superconductor sample while its imaginary component 
describes the rate with which the number of the Cooper pairs changes. 
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Figure 4: Real (upper at p « 0) and imaginary (lower at p ps 0) parts of the total radial 
functions (a) i?oo (p) and (b) i?oi (p) for several Aj where for both cases the thick solid lines 
are for Aj = 0.01, dash-dot-dotted lines - for A; = 1, and thin solid lines are for A, = 2 
while dotted lines are for (a) A» = 0.2 and (b) A, = 0.1, dashed lines - for (a) A, = 0.4 
and (b) Aj = 0.3, and dash-dotted lines - fojg(a) A, = 0.6 and (b) Aj = 0.5. 



3. Results and discussion 



In this section the results of the calculations based on the theory devel- 
oped in Section [2] are presented and their mathematical and physical inter- 
pretations are given. 

Fig. |2] shows Eqq and r 00 as functions of the pure imaginary de Gennes 
distance zA, and magnetic intensity B. It is seen that the real part of the 
energy as a function of A, at the fixed field decreases with the distance 
growing. As the field increases, the energy dependence on the extrapolation 
length gets smaller, and at the large field, as was deduced earlier, E^ m tend 
to the Landau levels, Eq. ( 129"|) . when the energy is A- independent and real, 
T = 0. However, at the zero and small fields, the imaginary part of the 
total energy exhibits a pronounced extremum on the Aj-dependence with 
it absolute value increasing at the decreasing field. Similar dependencies 
demonstrate also the states with other quantum numbers n and m. 

Before considering the magnetic field influence, let us discuss first the 
field-free case. In Fig. [3] complex energies E^ m are drawn as functions of 
the imaginary length z'Aj for several n and m. Panel (a) shows that for 
each state (n, m), as Aj changes from zero, the real part Enm monotonically 
decreases from {i\ m \ n /^) 2 at the pure Dirichlet case, A = 0, to (j| m | n /7r) 2 
characteristic for the Neumann requirement, |A| = oo. Here, j\ m \ n and j[ m | n 
are t he nth roots of the equations J\ m \{x) = and J!,(x) = 0, respectively 



131] . The panel inset makes a comparative analysis between the ground- 
state energy dependence on the positive Robin parameter A r and the real 
part of the transverse energy dependence on the purely imaginary de Gennes 
distance z'Aj of the same magnitude. It is seen that in either case the energies 
monotonically decrease to zero with |A| growing; however, in the former 
case the derivative dEoo/dA r is a monotonic A r function too while for the 
imaginary distances the derivative dE^/dAi has an extremum on the Aj 

(r) 

axis. Note also the faster approach to zero by i% for the large magnitudes 
of|A|. 

Panel (b) of Fig. [3] shows that the T nm — Aj dependence demonstrates 
highly nonmonotonic behavior with characteristic resonance whose sharpness 
and magnitude increase with n and \m\. This prominent feature can be 

understood as follows. An addition theorem and asymptotic properties of 

the Bessel functions [131| applied to ( 130]) lead to the following expressions 
for r n | m | in the limit of the small, |Aj| <C 1, and large, |Aj| 3> 1, magnitudes 
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of the purely imaginary de Gennes distance: 



r„H = — <; g^ngUJ j_ iai^i ( 46 ) 

5 nn / is a a Kronnecker symbol. Note that the coefficients at both Aj and 1/Aj 
are positive; accordingly, as the imaginary part of the extrapolation length 
moves away from the two extreme points Aj = and 1/Aj = 0, the magnitude 
of T nm grows. Matching these two asymptotics in the intermediate regime, 
Aj < 1, leads to the pronounced resonance with its maximum increasing 
with \m\ and n. Mathematically, this foll ows from the fact that j\ m \ n is an 
increasing function of both \m\ and n jl31 |. This increase is easily explained 



from the physical point of view too; namely, at the fixed principal number n 
(azimuthal number \m\) the states with the bigger \m\ (n) are located further 
from the origin; accordingly, they are influenced stronger by the change in 
the boundary conditions. In our model, such a resonance means that by a 
simple change of the applied gate voltage defining the field £ one can control 
an attenuation of the supercurrent over a wide range, as it follows from ( 14 ip 
and (|43|). 

Fig. H] depicts wave functions evolution as the purely imaginary de Gennes 
distance sweeps its positive axis. Nonzero imaginary part induces complex 
radial wave function what is a manifestation of the transformation of the 
true bound state with the infinite lifetime into the quasibound state with 
finite T nm . It is seen that the imaginary component of the wave function 
grows in magnitude for the small Aj, reaches its absolute maximal value 
approximately at the same de Gennes distance when the maximum of T nm 
is achieved and at the further growth of the extrapolation length gradually 
turns to zero again. At the same time, the real components transform from 
the Dirichlet dependence at Aj = to the Neumann case at the infinite 
de Gennes distance; in particular, the radial function Roo(p) approaches a 
constant value v2 for the large |A|. 

Current density patterns corresponding to the functions from Fig. H] are 
plotted in Figs. [5] and EJ They clearly show the emergence and the increase 
of the normal component of the supercurrent as the imaginary part of the 
extrapolation length reaches the resonance value and its gradual decrease 
to zero upon further increase of Aj. In the absence of the magnetic field, 
the state with m = for the two limiting cases |A| = and 1/|A| = 
does not show any transverse currents at all while for the finite Aj the radial 
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Figure 5: Current densities of the state with n = m = and (a) Aj = 0.01, (b) A,; = 0.2, 
(c) Aj — 0.4, (d) A, = 0.6, (e) Aj = 1, and (f) Aj = 2. Longer arrows denote larger 
currents. 
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component emerges, as Fig. demonstrates. For the nonzero azimuthal 
numbers the pure rotational whirling with j p = for the Dirichlet case is 
distorted at Aj ^ by the appearance of the radial motion which pushes 
carriers out of the sample (Fig. [6]). Similar to the m = case, this radial part 
reaches maximum approximately at the resonance extrapolation length and 
then with Aj growing gradually fades to zero transforming the picture into the 
currents in the Neumann disk. These current density patterns supplement 
theoretical discussion of the previous section about the carriers flow through 
the superconductor boundary. 

So far, we considered the case of the zero real part of the de Gennes 
distance. Fig. [7] justifies such an approach as its panel (b) shows that the 
finite positive A r strongly quenches the T nm resonance, and at A r > 5 the 
imaginary part of the transverse energy almost vanishes and is independent 
of A,. For the finite A r and small magnitudes of Aj, in the way, similar to 
f )46|) . one gets for the halfwidth T n \ m \: 



n\m\ 



7T 



(1 4" A r ) t/| m |(A| m | n ) + 2 A| m | n A r </| m |_| (A| m | n ) ^| m |_(_^ (A| m | n ) 



Aj 



where A| m | n is nth root of equation 



(47) 



I Ail <C 1, 



X 



0. 



For the vanishing real parts, A r — > 0, the solution Ai m | n tends to j\ m \ n , and 
fj47|) transforms into the corresponding limit of (|46p . In the opposite regime 
of |A r | — > oo the halfwidth approaches zero with the root Ai m i n changing to 
f\m\n- The expression for r n | m | in the limit |Aj| ^> 1 depends on the relation 
between A r and Aj and, in the scale of Fig. [?l basically coincides with the 
lower line of (|4*6|) . This is seen from Fig. [3(b) where for the large Aj the 
imaginary part of the energy only slightly depends on A r . Summing up 
this discussion, one concludes that the real background of the extrapolation 
length suppresses the influence of its imaginary component and restores the 
lossless particle flow without its significant changes in the external electric 
field. Accordingly, for the best experimental observation of the supercurrent 
attenuation or amplification, one needs to use ferromagnet as a surrounding 
material when, as stated in the Introduction, the real part of the de Gennes 
distance is very close to zero. The effect survives if the superconductor 
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;ure 6: The same as in Fig. [5]but for n = 0, m = 1 and (a) Aj = 0.01, (b) Aj = 0.1, (c) 
= 0.3, (d) A< = 0.5, (e) A, = 1, and (f) A. ( = 2. 
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borders the normal metal (0 < A r < oo) even though its magnitude gets 
decreased, and completely vanishes for the boundary with the insulator or the 
vacuum, 1/A r = 0, what is quite natural since in these ambient environments 
the Cooper pairs can not exist at all and, accordingly, they can not enter or 
leave the superconductor. 

As expected, the real part of the energy monotonically decreases to its 
Neumann counterpart if either part of the complex extrapolation length tends 
to the large values as is shown in panel (a). For the negative real parts of 
the Robin parameter the picture of the resonance quenching is basically the 
same except of the very vicinity of the negative zero when the energy takes 
unrestrictedly large negative values upon the extrapolation length approach 



to zero from the left [37|, |45| . In this case the increasing imaginary part pro- 



motes the negative energy Eq V J upwards to its zero Neumann case; however, 
at the very small negative A r , that energy could be occupied by the level 
with n — 1. Repulsive interaction between these two states leads to the an- 

(r) (r) 

ticrossing between i% and E[ on the Aj axis with its sharpness dependent 
on the smallness of the negative A r and results in the repel of the higher 
lying state from the lowest Neumann position. To make our exposition more 
compact, we do not present this transformation here. 

It is instructive to compare Figs. |2] and [71 Panels (b) show quite a similar- 
ity, at least qualitative, while the real parts of the transverse energy exhibit 
completely different behavior. Thus, in order to observe the most profound 
supercurrent alteration, one needs to use materials with A = and, to check 
the quenching of the resonance, it is necessary to turn the external magnetic 
field on. 

The last paragraph naturally returns us to the study of the combined 
influence of the magnetic field and complex extrapolation length. We have 
already seen in Fig. MJd) that the increasing magnetic field washes out the 
resonance on the T — A, axis since it pushes the Cooper pairs closer to the 
wire axis and thus restricts their interaction with the disturbing borders of 
the wire. For the very large magnetic fields all the states (n, m) approach 
from below the corresponding Landau levels, Eq. ( 129]) . what means their 
independence on the imaginary component Aj. It is known that for the 
real extrapolation length similar asymptotics of the Landau levels crucially 
depends on A r ; namely, it was calculated that for the Dirichlet confinement 
the state with m = at any magnetic intensity remains the lowest one of 
the whole family of the levels with the same principal quantum number n 
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1221 . Il23| . Situation is very different for the Neumann case, 1/A r = 0, where 
the levels with the adjacent negative azimut hal quantum n umbers —\m\ and 
— \m\ — 1 cross with the field increasing 71, 1261 . 129 . 13dl | thus causing the 
ground state to change its vorticity. These theoretically predicted crossings 
were ind eed observed experimentally for the mesoscopic aluminium disks 



14lHl43j . The crossings are characteristic for any nonzero extrapolation 
length with their location shifted to the higher magnetic fields for the smaller 
de Gennes distance until at A = they disappear and the state (0, 0) remains 



the ground state for any B [7l|]. Fig. [8] showing the dependence of the real 
part of the transverse energy on the magnetic field for the purely imaginary 
distances Aj = 0.3 and A, = 0.5 confirms this property for the complex 
extrapolation lengths too. These changes by the ground state of its azimuthal 
number are responsible for the jumps of the magnetization. 

Magnetizations M z for the states (0, 0) and (0, —1) are shown in Fig. [9]as 
functions of the magnetic field and imaginary de Gennes distance iAj. The 
magnetic moments from its zero-field value — m [see Eq. fl37|) for B = 0], 
decrease for the large B to the magnetization of the free charged particle in 
the uniform magnetic field M z = — (2n+ \m\ +m + 1). For the Dirichlet disk 
this approach is a monotonic function of the intensity while for the Neumann 
case the magnetization reaches minimum as a function of the field after which 
it monotonically increases to the saturation value. Extremum location on the 
B axis increases with the azimuthal number \m\. A transformation between 
the two cases with Aj changing is clearly seen in the figure. Contrary to the 
longitudinal supercurrent features, the dependencies in Figs. |8] and [9] are very 
similar to the case of the real de Gennes lengths. Accordingly, magnetization 
measurements are hardly suitable for the detection of the complex boundary 
conditions. 



4. Concluding remarks 

The analysis of the model of the superconducting 2D disk or infinitely long 
3D cylinder with the nonzero imaginary part of the Robin parameter con- 
firmed the general quantum mechanical rule that the complex boundary con- 
dition leads to the complex eigenenergies of the corresponding Shcrodinger 
equation and induces currents through the superconductor interface into the 
surrounding medium. For the case of the infinite wire these transverse leaking 
currents alter the longitudinal flow of the Cooper pairs with the resonance 
dependence on the imaginary part of the de Gennes distance. 
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Figure 8: Real parts Eq^ of the total transverse energies EQ m as functions of the magnetic 
field B for Aj = 0.5 (upper panel) and Aj = 0.3 (lower panel). Numbers near the curves 
denote corresponding azimuthal quantum numbers m. 
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Figure 9: Magnetizations M z of the states (0, 0) (upper panel) and (0,-1) (lower panel) as 
functions of magnetic field B and purely imaginary de Gennes distance iAj. Note different 
M z scales for the two panels. 
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As was stated in the Introduction, the present treatment for the zero 
magnetic field can be directly applied to the study of the dynamics in the 
acoustical, impedance and ferrite-filled electromagnetic waveguides as well 
as the processes in the superconducting structures close to the transition be- 
tween the normal and superconducting states. The superconductors in the 
magnetic field are also described by the theory developed above. As was 
mentioned in Sec. 121 for the thin superconducting disks the results obtained 
are sup porte d by the nonlinear GL theory as well. Moreover, numerical sim- 
ulation 144j | of the aluminum disks with the Neumann requirement revealed 
that the predictions o f the linearized GL equation are in a better agreement 



with the experiment [14311 tha n the full GL treatment. Furthermore, theo- 



retical linear GL analysis [144] correctly captures experimental findings [143 
down to temperatures T = 0.4 K lying well below corresponding T c = 1.19 K 
for aluminum. Accordingly, there is a strong confidence that the resonances 
discovered in the present investigation, could be experimentally detected in 
the superconducting guiding structures with the applied external electric field 
tuning the total complex extrapolation length A. 

The de Gennes boundary condition uses only one parameter A for the 
description of the interaction between the superconductor and surrounding 
medium and ignores the processes occurring outside the superconducting 
sample. However, the Cooper pairs can penetrate into the adjacent normal 
metal and, obviously, into the neighboring superconductor with the higher 
critical temperature. This can be especially important in our case of the 
complex boundary conditions when these charge carriers flowing outside the 
actual superconductor can contribute to the total longitudinal supercurrent 
and magnetization. Modified GL theory taking into account nonz ero o rder 



parameter in the bo rdering normal metal was recently proposed |l45j and 



developed |146Hl49| . Instead of the single parameter A, it utilizes the two 
ones: the carrier mass rh (here we operate in regular, dimensional units and 
use tilde for the values of the modified GL theory) which is different in the 
normal metal with its value of rh n and in the superconductor (m s ), and the 
dimensionless parameter a = \(m n a n ) / (rh s a s )\ where the GL parameter —a 
has different signs inside and outside the superconductor, ot s a n < 0. As a 
result, it requires on all the interfaces a continuity of the vector potential A, 
order parameter ^, and values of n4- (iti^ — qA^j ^ and i|vx x n 

with /i being the global magnetic permeability. Taking into account the 
currents flowing outside the superconductor can alter the results presented 
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here. However, this modification should not smear out the resonances at 
all, especially for the case of the bordering ferromagnets which effectively 
destroy superconductivity by breaking the Cooper pairs. In a sense, the 
model of the complex extrapolation length is similar to the modified GL 
theory since both of them use the two parameters for the description of the 
superconductor /other medium interface (in our case these are the real and 
imaginary components of the de Gennes distance) and allow the current flow 
across the border. 

For considered here system of the quantum disk there is only one surface 
through which the energy can flow into or out from the structure. In turn, 
for the quantum ring with the interfaces of the smaller and larger radii there 
are two such channels. Intuitively, one can expect that for some correlation 
between the inner and outer complex de Gennes distances with the opposite 
signs of their imaginary parts the energy can become real again. Apparently, 
this happens when the flux entering through one of the surfaces is equal 
in magnitude to the outgoing flow through the other boundary. For the 
rectangular geometry it is elementary to show that it happens when the de 
Gennes lengths on the opposite walls are complex conjugate (see E g. (1 0) in 



Ref. |45j). However, for the annulus due to the centrifugal forces [150j this 
relation does not seem to be so trivial. This case requires separate special 
consideration. 
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